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^ ! Abstract 

In a recently proposed multi-Higgs extension of the standard model in which discrete symmetries, 
■ A4 and Z3 are imposed we show that, after accommodating the fermion masses and the mixing 

CN ■ matrices in the charged currents, the mixing matrices in the neutral currents induced by neutral 
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scalars are numerically obtained. However, the flavor changing neutral currents are under control 
mainly by mixing and/or mass suppressions in the neutral scalar sector. 
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I. INTRODUCTION 



One of the main motivations to go beyond the standard model is to have some hints 
about the flavor problem, that is, to understand the pattern of fermion masses and mixing. 
Most electroweak models have mass matrices of the form M^p = J2i(Xi)ai3{^'i) y where 
the FjS are, for Dirac fermions, arbitrary complex dimensionless 3x3 matrices, and ($°) 
denotes the vacuum expectation values (VEVs) of the neutral scalar fields in the model. 
For Majorana fermions, the FjS are complex symmetric matrices. The mixing matrix and 
the mass pattern in each charge sector depend on the structure of the respective FjS. It is 
well known that explicit and predictive forms of these matrices can be obtained by imposing 
flavor symmetries. In particular, most of the ansatze for the fermion mass matrices are of 
the form |l| 
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and other textures which have in common the fact that the zeros appear in symmetric 
entries. Recently, it was proposed an extension of the electroweak model with A^and 
discrete symmetries, in which the quark mass matrices have the following texture 



^ Aq Dq^ 
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where the zeros occur in non-symmetrical entries. 
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II. QUARK MASSES AND MIXING IN MODEL WITH A4 SYMMETRY 

We consider a model with 6*57;/ Z3 (g) Z3 (g) Z3 symmetry^here Gsm is the standard 
model (SM) gauge symmetry, SU{3)c ® SU{2)l ® U{1)y- In [2] quarks and leptons were 
considered briefly. Here we show a more detailed study in both sectors. Let us begin 
by considering quarks. The model needs twelve SU{2) doublets, four triplets of the ^4 
symmetry: H = {Hi, H2, H^), and $ = {^1,^2,^3) for the masses of the M-quarks; H' = 
{H[,H2,H^) and $' = {^[,^'2,^'^) for the masses of the d-quarks. We also introduce a 
complex scalar singlet ( [2]- 
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Fermion weak eigenstates transform under {A^, Z3, Z'^, Z'^) as: 

Ql = {QiL, Q2L, Qsl) ~ (3, c<;^ 1, 1); Ur = {u2r, usr, uir) ~ (3, 1, 1, w); 

Vr = {d2R,d2,R,diR) {Z,l,uj,uj), (3) 

in which oj = e^'^*/^, and Qn = {di Ui)]^, i = 1,2, 3. All fermion fields in are symmetry 
eigenstates. Notice that, assuming the usual ordering of QiL, the A4 symmetry will allow, 
after its breaking, to distinguish among the right-handed quark components. The scalar 
fields transform under the same symmetries as: 



H ~ (3, u, 1, u), H' ^ (3, w^), $ ~ (3, u, u, u), 

<l>'~(3,u;^l,u;2), C~(l,l,c^,l). (4) 

The notation in the scalar sector is slightly different from that in ^ . 

With these fermion and scalar fields, and the discrete symmetries above, we have the 
leading Yukawa interactions in the quark sector 



hJQ^HU + h'^^[Q^^R]UR 



+ 



-H.c{5) 



where (p = ^4'*, with e = ia2 (o"2 being the usual Pauli matrix); [Xy]^^^ means the appropri- 
ate product, and denotes the singlet state, see Ref. and references therein. Instead 
of the A and B triplet representations we can use the symmetric, s, and antisymmetric, a, 
ones. In fact, A = s + a and B = s — a. We choose to have some fields which products are in 
A (or s), and others in B (or a). We assume this, firstly, because in global symmetries it is 
not mandatory to use all the product representation allowed by the symmetry and, secondly, 
because that choice may be explained by an underlying still unknown dynamics. Notice that 
the symmetry imposes only two dimensionless Yukawa couplings in each charge sector. 
We assume that V/A ^ where V denotes any VEV of the model and A is an energy 
scale characterizing an unknown physics. Notice that the model has flavor changing neutral 
currents (FCNC) since several Higgs doublets contribute to the masses of a given charge 
sector {4]. 

The renormalizable interactions in (|5]) induce diagonal interactions in the weak basis 
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which exphcitly read 



+ [Q2LH3d2R + Q^LH'id^R + Q^LH2dm] + H.c. 



(6) 



All scalar doublets are of the form [xf On the other hand, the non-renormalizable 

interactions ([5]) are written explicitly as 



— r^' 



= h'^j[Q3L'^2U2R + QiL$3M3fl + Qai'^l^li?] 
+ KjiQ3L^2d2R + QlL^sdsR + Q2L^[diR\ + H.C. 



(7) 



The mass matrices obtained from Eqs. ([6]) and ([7]) are 



V2 











^^1 



+ H.C. 



(8) 



for the u-type quarks, and 



V 











+ H.C., 



(9) 



for the d-type quarks. The texture of these matrices are different from those of other multi- 
Higgs models like the private Higgs in j^, and we stress that they are a consequence of 
the Ai symmetry and, mainly, of the choice of the representation A or B. Above we have 



h' 

defined a,=^^ and a 



h' 

■jf-^-i for the 2/3 and —1/3 charged quarks, respectively, and 



we have denoted = fj, (0°) = f^., (/if) = f •, = f^., and {Q = f(. For the sake of 
simplicity, all parameters in Eqs. (|8]) and iQ have been considered real. Notice that since 
each charged sector has its private VEVs, these matrices are independent from each other. 

Here we will show, numerically, that with the mass matrices ([8]) and (|9]) it is possible to 
accommodate the observed masses and the mixing matrices in the quark sector. These mass 
matrices are diagonalized by bi-unitary transformations, Vj^ M^V^^^ = and V/^M^V^^ = 
M^, respectively, with = diag{mu,mc,mt) and = diag{md,ms.,raij). The change 
of basis is qiL{R) = {V^^-j)iaqaL{R), Qa denotes the quark mass eigenstates of the respective 
charge sector, = u,c,t for quarks with electric charge 2/3 and qa = d,s,b for for quarks 



with electric charge —1/3. We obtain the unitary V^'^ matrices, by solving the matrix 
equations: 

M.Mtvf t = dtagimlmlml) = (MJ^ V^MIM^V^^ = (MJ^, 
M.Mjvf^ = dtagimlmlm',) = (MJ^ V^MIM^vP = (MJ^ 
lfV« = l, VS^V2 = 1, Q = U,D, (10) 

using as input parameters those VEVs, and hg, appearing in M^,M^, which give the 
observed quark masses at an appropriate energy scale (the Z mass in the present case) and 
the Cabibbo-Kobayashi-Maskawa quark-mixing matrix, defined as V^^^^^ = Vi^V^^ . The 
latter one is such that, at least in the context of the SM, the mixing between the first two 
families is almost scale independent and for the third family Vet and Vub change at the level 
of 13-16% between rrit and 10^^ GeV Q]. Here, therefore the V^^j^/ ^^^^ considered scale 
independent. 

Since the mass matrices ([8]) and (Q are predictions of the model, they are valid at the 
energies at which all symmetries of the model are realized, i.e., at the electroweak scale. 
For this reason, we use the running quark masses at /i = Mz, taken from Ref. j^], for light 
quarks (in MeV): m„ = 1.27lg;42, ''^d = '^-QOtl'^g, mg = 55tll; heavy quarks (in 

GeV): = 0.619 ± 0.084, rub = 2.89 ± 0.09, rut = 171.7 ± 3.0. 

In order to obtain the values of the V^'^ matrix elements within an interval, we find two 
sets of values of the input parameters which give the quark masses and the CKM entries 
within the experimental errors. For instance, using: 1) /i^ = 0.1 and = 0.11 and (all 
VEVs are given in GeV) v[ = 28.9, v'^ = 0.44, v'^ = 0.03, v'^^ = 1.95, v'^^ = 0.03, and 
17^3 = 8.1, and 2) = 0.1 and = 0.2 and v[ = 29.8, v'^ = 0.54, v'^ = 0.03, v'^^ = 0.35, 
~ 0.01, and v'^^ = 2.9, we obtain the following values for the masses in the d-quark 
sector: = (2.70-2.97) MeV, m, = (48.95-54.4) MeV, and = (2.89-2.98) GeV. For 
the M-quark sector, we use 1) = 1.11, = 0.2, and the VEVs Vi = 153, V2 = 0.54, = 
0.001125, v^^ = 0.08875, v^^ = 0.03555, v^^ = 57.2355; 2) h^, = 1.11, = 0.13, and the 
VEVs vi = 153, V2 = 0.531, V3 = 0.00108, v^^ = 1.2048, v^^ = 0.3199, v^^ = 70.783, we 
obtain: m„ = (1.27 - 1.93) MeV, = (598 - 613) GeV, and mt = (170.305 - 170.137) 
GeV. In spite that some VEVs are small, this does not imply necessarily the existence of 



light scalars. See Ref. js] where the case of three Higgs scalar doublets with A4 symmetry 
were considered in details. 
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The matrices V^'j^ can be obtained by solving ffTOl) with values for the parameters con- 
sidered above. An example of numerical V^'^ matrices obtained is given by: 



V, 



u 



( 



D 



\ 



0.03290 ^ 0.28285 -(0.95902 ^ 0.99946) -(0.00246 ^ 0.01678) 
-(0.95733 ^ 0.99667) -(0.0330 ^ 0.28334) 0.05680 ^ 0.07457 
0.05923 ^ 0.07461 (0.3 2) x 10"^ 0.99721 -> 0.99824 

0.12895 ^ 0.25492 -(0.96693 ^ 0.99165) -(0.00251 ^ 0.00786) 
-(0.96647 ^ 0.99146) -(0.12896 ^ 0.25504) 0.0193 ^ 0.0298 
0.01946 ^ 0.03082 (2 ^ 9) x 10"^ 0.99952 0.99981 



Using these matrices and the definition V^^^^^ 



V^Vl^\ we obtain 



\ 



0.9748 - 0.9875 0.1571 -0.2230 0.0014- 0.0113 
0.1574- 0.2227 0.9739 -0.9868 0.0381 - 0.0438 
0.0051 - 0.0111 0.0394-0.0424 0.9990 - 0.9992 



(12) 



We compare this matrix with the global fit of the magnitude of the V^^^j^j elements given in 
Eq. (11.27) of PDG Q 

( 0.97428 ±0.00015 0.2253 ± 0.0007 0.00347+^-^^^}^ \ 



\ypdg 



0.00012 

0.2252 ± 0.0007 0.97345l[!:J^[!J^}^ 0.010^J^;[!|^^^ 



V 



0.00862l°:S 



0403+°-°°^^ 

U.U^UO_o.0007 



9991 c;2+° °°°°3° 



(13) 



Considering the elements of V'J'^^^^, we note that the V^j^j^j in Eq. f[T21) has all its entries 
within l-cr but the Vtd-, which is within 1.7-a. We recall that at present there are several 
discrepancies between experiments and the standard model at the tree level that are about 
3(7 standard deviations [1C|. Notice that we are not considering CP violation in our analysis. 
This is because in our framework the CP violation issue is more complicated than in the SM 
case since in our model there are more CP violation phases. Although we can perform the 
usual phase redefinition in the left-handed fields of the charged currents, in order to obtain 
the full, single-phase, Vckm niatrix, many CP violating phases will still be present in the 
Yukawa interaction Lagrangian. We find that this subject deserves a separate study. 
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In the same way we have obtained, from f fTOj) . the respective Vj^'^ matrices: 



V, 



R 



\ 



J 



V, 



D 



-(0.00091 0.00371) -(0.99999 ^ 1.) 0.00005 ^ 0.00027 
0.99999 ^ 1. 0.00091 -> 0.00372 -(0.0002 0.00026) 
0.0002 ^ 0.00026 0.00004 ^ 0.00027 1. 

( 0.00047 ^ 0.01661 -(0.99986 ^ 0.99997) 0.00006 ^ 0.00011 ^ 
-(0.99986 ^ 0.99997) -(0.00703 ^ 0.01661) 0.00035 ^ 0.00047 
0.00035 ^ 0.00047 0.00007 ^ 0.00011 1. 



(14) 



The matrices above will appear in flavor changing neutral currents in the Yukawa interac- 
tions. 

Notice that the mass matrices in Eqs. ([8]) and (Q have the following parameters (assuming 
all of them to be real): h^,a^,h^,a^ and twelve VEVs. It means 16 real free parameters 
to explain 12 mixing angles of four matrices: V^'^ (or one of V^'^ and V^j^)aTid V^'^. 
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12| in order 



However, if CP violation is allowed, it is necessary to chose a weak basis 
to eliminate some extra phases which may be difficult to measure. 

FCNC suppression. There are flavor changing neutral currents (FCNC) effects in both 
quark sectors. Here we will consider only some of these sort of effects which occurs at the tree 
level. For instance, from the renormalizable interactions of the doublets with the d-type 
quarks (in the weak basis), are given by h ^{d2LS2R'{''-i + d'udmVi + '^ilc?irV^2')- The other 
three doublets have interactions, which are suppressed by the scale A, will not be considered. 
On the other hand, the neutral scalar weak eigenstates ^p'^ are linear combinations of the 
neutral scalar mass eigenstates, i.e., ip'^ = J2n^anhn, a, n = 1,2,3. This is in fact 
a simplification, since the model has several doublets, but it may be assumed (for the 
sake of simplicity) that in each charge sector there are three doublets which give the more 
important effects. It means that the mixing matrices in the full scalar sector may be almost 
block diagonal, with each one related to a given fermion charge sector. 

The suppression of FCNC can be obtained at least with a reasonable fine tuning in the 
mixing parameters or by considering heavy enough the charged and neutral scalars. We 
consider some examples just for illustrating this point. The contributions of the neutral 
scalars to the AM^ put the strongest constraint on some of the parameters of the model. 
From dnj, the renormalizable Yukawa interactions between the d and s quarks and a given 
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neutral (pseudo) scalar, denoted by are given by 

- Cds = J^l^LKnSR + SLKdR]{hl + iAl) + H.c. 
n 

where we have assumed /i„ real and defined 



(15) 



2n 



^ -0.01f/3n + 0.16?72„ + 10-^f/i„, 

K = {vE):siyE)sdU,n + iyE)iiyE)uU,n + {yEMyE)MU2n 

^ O.Olt/sn - 0.01^72n + 10~^[/i„, 



(16) 



where we have used the values of the VE and VE matrix elements given in ffTTj) and ([HD, 
respectively. 

Thus, the interactions in f lT^ can be rewritten as 



= ^{[{K^ + K':){ds) + (ir„ - ir::)(j75s)](/^° + ^O] 
+[(ir: + K:)(^rf) - {k: - K^s^.d^hi + ^o*}. 



(17) 



and the effective Hamiltonian contributing to K ^ K transition is given by 

■H^ff=%calars = ^^^^n + Kf{sdf + {K - K)\s^,df] (18) 

n 

In the vacuum insertion approximation {K^\{sd)'y_ji\K^) = 2MKfK/S. Using (up to 
some phases) Q 



(ir°|(sc/)2|K°) = 
{K^\{sj,dY\K') 



PkMk 



12 



Ml 



K 



12 



1 - 11. 



Ml 



{rus + may 



we obtain the following extra contributions to AM^ due to the neutral scalar 



(19) 



AM, 



K scalars 



2Re{K'\ntffXcaiars\K') = Re Ott^x/. 



2 



(20) 



where 



sd 



/)2 



n 



1 - 



Ml 



7.52(f/3*„ + f/3„)t/2n + 0.75f/2 



(m„/100GeV)2 



-2 

2" 10"^ GeV"^ 



+ {K - K 



I \2 



1 - 11- 



Ml 



(21) 
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where we have used only the main contributions in ( !T6|) . There are also similar contributions 
induced by the pseudoscalar, A^. For illustrative purposes we showed above only the scalar 
contributions. 

Notice that, in the present context, in flTB]) only the matrix elements Uan are not known 
yet and that, independently of the mass eigenstates Uin will not be constrained by 
processes like AMk- In order to be consistent with data Re (^^ must be smaller than the 
contribution of the SM: i.e., ReCf/ = GlmlRe[{VcKM)*diycKM)cs]VlQTr^ ~ IQ-i^GeV"^ 
(we have used only the dominant contribution of the c quark and g{mc/Mw) = !)• By 
imposing that Re J2nQd < ^^Csf^y implies, from f ET]) 

IdReUsn ■ ReU2n + 0.75 [(Ref/2„)' - {lmU2n?)] 



< 10'". (22) 



(m„/100GeV]2 

Even if for a given n, rrin = 100 GeV, there are two different ways for each term in (22) 
to satisfy the constraint: a) all parameters involved are of the order of 10~^, i.e. KeU^n ~ 
Re?72n ~ Inif/2n ~ 10~^, or b) there is a fine tuning among the parameters. For instance, 
ReUsn = 0.01, Re [/2n = 0.01,Imf/2n = 0.045826742. For heavier neutral scalars the Uan 
matrix elements may be greater, limited only by the unitarity of the matrix. 

Similar constraints come from AMb and AMb^ data and other AB = 1 weak 
processes, like Bg^d — > A*^/^^ and -^'''|-| ~^ K~^fi~^fi~ decays. Experimentally, 
T{B^ — 7- fi^ fi')/T total < 5.2 X 10^^ [9]. In this case, the quark interactions in- 
volved are h,h [iV£rsM)ss^f + iV£UV£),s^f + iVf^YMU^f] sn + H.c, and we 
recall that, in this model, leptons have their own scalar sector with the Yukawa inter- 
actions gi[[LH']AlFt]i- Here, H denotes the triplet of A4 formed by three scalar dou- 
blets (see below). For muons it means giJIi^fi^H^, and the muon mass is given by 
= gi{h^ = giv^. In the scalar potential the scalar triplet H mixes with the scalar 
triplets related to the quark sector only by quartic terms like Ai H^H [if'^if']^ = 



Ai 



A 

This implies terms like Aif2'03 in the mass 



matrix in the neutral scalar sector and, in the weak basis a mass insertion implies that 
the propagator becomes Aifg-Os/m^, where m denotes a typical value for the neutral scalar 
masses. As we said before, in that basis, the Higgs scalars coupled to leptons are different 
from those coupled to quarks. It means that in semi-leptonic decays a detailed analysis 
may be done only if we consider at least six doublets of scalar Higgs bosons. Here, just for 



9 



illustration consider b — )■ s/i/i decay. Working in the weak basis we obtain, 



oc 



2 2/2 

^ = 1.155 X 10^A?/(m/GeV)^ < 10"^ (23) 



r(fe CZ/^/i) GpTTl^ 

and we have used the mixing matrix elements (yl')db and (y^)ds from ( !TT|) and ( fT4l) . 
respectively. All parameters in ( l23l) but Ai and m are already known. It implies 
A^/(m/GeV)® < 10^^^ which is valid for a wide range of the parameters. For instance, for 
m > 10 GeV and Xf ~ 10"'^ the above condition is satisfied. 

The case of FCNCs in the w-quark sector involves the matrices Vj^ and Vj^ . The ana- 
lysis of the FCNC in this model is similar to that of Ref. [l6| in which two doublets were 
considered. 

The lepton sector. In the lepton sector the discrete symmetries are the same as in the 
quark sector. Leptons transform under A4 ® (Z^)^ as follows: 

L = (Le,L^,L^) ~ (3,w,a;^ 1); Ir = {hr,tr, cr) (3,1,1, u), (24) 

and the leptophilic scalars transform under those symmetries as 



H" = [HlHlH'^ ~ (3,u;,u;^l), H = {HuH2.H^) ~ {Z,uj,uj\uj^), 

$" ^ ($;', $'2', $3) ~ (3, 1, 1, 1), r = (Ti, T2, Ts) ~ (3, uj\ 00, 1). (25) 

The Lagrangian of the lepton sector is given by (see the notation in j^): 
£ = (^g^[LHU + jL[LHU\C\' + ---yR 

+ I (/. [(UeH")UiLe<^")]B + ^[L^ea ■ f]B[L<l>'%x + ..•)+ H.c, (26) 

where gi, g[, fv and fl are dimensionless Yukawa couplings and A an energy scale which may 
be, or not, the same as that in the quark sector. The mass matrix for the charged leptons 
is almost diagonal: Mi = gidiag{v2,v^,Vi) + 0{g[v'^Vil Is?'), where {Hi) = Vi. In this case 
the renormalizable interactions are dominant, i.e., gfii ~ m^, giv^ ~ e giV2 — me. 
Neglecting the contributions proportional to g[, the charged lepton mass matrix is diagonal, 
thus the values of the VEVs are easily obtained: Just for illustrating, by using the central 
value for the lepton masses (in MeV) at the Z pole scale: me ~ 0.486, 102.718 and 

rrir ^ 1746.24, hence V2 ~ rrie/gi, h ~ rn^/guVi ~ rrir/gi. 
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On the other hand, the mass matrix for the neutrinos is 



't>3 



\ 



"2 X "3 'P2 



03 "3 n X 



A 



1 5 



(27) 



in which 6i is given by 



5i 



fu 



fu Afif</,i 



(28) 



where (x) = v^, e (A°) 



Since the mass matrix for the charged leptons is, for practical purposes, diagonal, which 
imphes ?7| ^ = 1, the mixing matrix PMNS, is defined as Vpmns = = U, i.e., is obtained 
directly from the diagonalization of the neutrino mass matrix. 

In the neutrino sector the Vpmns niatrix compatible with experimental data is, at the 
3a level llTl], 

^0.77- 0.86 0.50 -0.63 0.00 - 0.22 
0.22 - 0.56 0.44 -0.73 0.57- 0.80 
0.21 - 0.55 0.40 -0.71 0.59 - 0.82 
We use a different strategy from that used with quarks. Imposing that the matrix U satisfies 
the equation U^M^U = Mu = diag(mi, m2, ms), where Mi, is given by the matrix in 
( 127|) . we found numerical values for all the parameters in this matrix. Next, we found 
the square mass differences. For example, using the lower limit of the entries in ( I29l) . and 



3a 



\ 



J 



(29) 



the values of A = 1 TeV, = 2 GeV and vt^ 



2vt3 = 1-25 GeV and (in eV) 



4,f2 = W\v3 = 0.04, w,^, = 40450, = 48000, t;^3 = 34850, we obtain (in eV) 



v 



0.03485 ~ 
~ 0.0337083 ~0 
~ 0.06 



(30) 



where ~ means entries smaller than 10 eV. From fl30|) . the mass squared differences (in 



eV ) obtained are 



Am2i 



8 X 10"^ Am^i = 2.5 x 10~^ (31) 

which agree within la with the values at the /i = Z given in Ref. The neutrino sector in 
the model differs from that of the private Higgs of [18^ but, both models can accommodate a 
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nonzero ^13 angle. Notice that if we assume that L{H") = L{T) = —2, and the other scalars 
having L = 0, the interactions in f l26|) conserve the lepton number. However, there is no 
Majoron since the ^4 and Z3 symmetries allow terms in the potential which break explicitly 
the lepton number. For instance, [H"^^]i[H^ H]i, [He^']i[H'eH"]i, H"]i[HeH']i. As in 
the quark sector, there are also FCNCs in the lepton sector but they will be considered 
elsewhere. 

Conclusion. Motivated by the different mass scales in the quark and lepton sector we 
propose a model in which each charge sector has its own Higgs scalars which acquires ap- 
propriate VEVs. In this case, in order to induce the respective fermion masses, the Yukawa 
couplings do not need to have a large hierarchy among them, i.e., all of them may be of 
the same order of magnitude The hierarchy is translated to the values of the VEVs but 
these could be, in principle, explained by the minimization of the scalar potential. In fact, 
it was shown in jsj that this is indeed the case, for three scalar doublets. 

Matrices like those in ([T]) are written in terms of dimensionless parameters. This is 
because in the context of the SM, where these ansatze have been usually studied, the mass 
scale is already determined by the VEV f^^^ = ((7^/2-^/2)^/^ but in multi-Higgs models 



only \/^~vf which satisfies this constrained. This is the case of the model in Ref. [2| that 
we are considering here, and all VEVs are considered parameters to be fixed only by the 
fermion masses in each charge sector. We recall that we are not considering physical phases, 
but Vj^'^ may have the number of phases that are allowed for an arbitrary unitary matrix. 
We allow that in the V^^j^,^ matrix only one physical CP violating phase to survive after 
redefining the quark fields. 

Another concern is if the suppression of the FCNC effect at low energies is stable up 



to higher energies. The answer is yes, at_ 
model 



order 



east for the case of a general two Higgs doublet 



. Moreover, QCD corrections [20| have also to be considered at the next-to-leading 



22| . This has been done in the 2HDM in the context of natural flavor conservation 
and minimal flavour violation in Ref. 23| . A similar analysis in the context of models with at 
least three scalar doublets with or without extra symmetries, say ^4, 83, will be considered 
elsewhere. 
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